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We consider several classes of symmetries of the Dirac Hamiltonian in 3+1 dimensions, with 
axially-deformed scalar and vector potentials. The symmetries include the known pseudospin and 
spin limits and additional symmetries which occur when the potentials depend on different variables. 
Supersymmetries are observed within each class and the corresponding charges are identified. 
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The Dirac equation plays a key role in microscopic 
descriptions of many-fermion systems, employing covari- 
ant density functional theory and the relativistic mean 
field approach. The relevant mean-field potentials are of 
Coulomb vector type in atoms, and a mixture of Lorentz 
vector and scalar potentials in nuclei and hadrons [lj. 
Recently, symmetries of Dirac Hamiltonians with such 
mixed Lorentz structure have been shown to be rele- 
vant for explaining the observed degeneracies of certain 
shell-model orbitals in nuclei ( "pseudospin doublets" ) Q , 
and the absence of quark spin-orbit splitting ( "spin dou- 
blets") as observed in heavy-light quark mesons. 
Supersymmetric patterns have been identified in spe- 
cific limits of such spherical potentials 0, 0]. In the 
present Letter we further explore classes of symmetries 
and supersymmetries when these potentials are axially- 
deformed. Such a study is significant in view of the fact 
that mean-field Hamiltonians often break the rotational 
symmetry. Cylindrical geometries are relevant to a num- 
ber of problems, including electron channeling in crystals, 
structure of axially-deformed nuclei and quark confine- 
ment in spheroidal flux-tubes. 

The Dirac Hamiltonian, H, for a fermion of mass M 
moving in external scalar, Vs, and vector, Vy, potentials 
is given by H = a ■ p + p(M + Vs) + V v |{. When the 
potentials are axially-symmetric, i.e., independent of the 
azimuthal angle (p, Vg,v = Vs,v(p,z), P = \/x 2 + y 2 , 
then the z-component of the angular momentum opera- 
tor, J z , commutes with H and its half-integer eigenval- 
ues il are used to label the Dirac wave functions ^ = 



*l 2 , if+e- 
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Here 



g ± = g ± (p,z) and f ± = f ± {p,z) are the radial wave 
functions of the upper and lower components, respec- 
tively. Henceforth, such a wave function will be denoted 
by : {<7 + > 9~i f + i / }• The potentials enter the Dirac 
equation through the combinations 

A(p,z) = E + M + V s {p,z)-V v (p,z) (la) 
B(p, z) = E — M — V s (p, z) - V v (p, z) . (lb) 



bound Dirac states with \E\ < M and normalizable wave 
functions in potentials satisfying pVs(p, z), pVy(p, z) — > 
for p — > and Vg(p, z), Vy(p, z) — > for p — > oo or 
z —* ±oo. The boundary conditions imply that the ra- 
dial wave functions fall off exponentially for large dis- 
tances and behave as a power law for p — ► 0. Further- 
more, for z = and p — > oo, f~/g + oc (M — E) > 
and fir// + oc (M + E) > 0, while for z = and p 0, 
/~ /g + oc B(0)p and g~ / f + oc — A(0)p. These properties 
have important implications for the structure of radial 
nodes. In particular, it follows that for potentials with 
the indicated asymptotic behaviour and A(0), B(0) > 0, 
as encountered in nuclei, a necessary condition for a node- 
less bound eigenstate of a Dirac Hamiltonian is 
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The Dirac equation, — E^, leads to a set of four 
coupled partial differential equations involving the radial 
wave functions. Their solutions are greatly simplified in 
the presence of symmetries. We now discuss four classes 
of relativistic symmetries and possible supersymmetries 
within each class. 

The symmetry of class I, referred to as pseudospin 
symmetry, occurs when the sum of the scalar and vec- 
tor potentials is a constant, Vs(p, z) + Vv{p,z) = Ao. 

The symmetry generators, Si, commute with the Dirac 
Hamiltonian and span an SU(2) algebra 0,11 
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Here Si = Ci/2 are the usual spin operators, defined in 
terms of Pauli matrices. The Dirac eigenf unctions in the 
pseudospin limit satisfy 



fx = ±1/2 



(4) 



and form degenerate SU (2) doublets. Their wave func- 
tions have been shown to be of the form 19] 



For each solution with > 0, there is a degenerate time- 
reversed solution with — O < 0, hence, we confine the dis- 
cussion to solutions with Q > 0. Of particular interest are 



n 1= A-i/2 ' xy ' - h ' J > ' 
n 2 =A+i/2 ■ \y> y > Ji U J > 



(5a) 
(5b) 
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where A = ft — /i > is the eigenvalue of J z — S z . The 
relativistic pseudospin symmetry has been tested in nu- 
merous realistic mean field calculations of nuclei and were 
found to be obeyed to a good approximation, especially 
for doublets near the Fermi surface [3, [lfj. The domi- 
nant upper components of the states in Eq. |5|), involv- 
ing g + and g~ , correspond to non-relativistic pseudospin 
doublets with asymptotic (Nilsson) quantum numbers 
[N, n 3) A] ft = A + 1/2 and [N, n 3) A + 2] ft = A + 3/2, 
respectively. The doublet is expressed in terms of the 
pseudo-orbital angular momentum projection, A = A + l, 
which is added to the pseudospin projection, p = ±1/2, 
to form doublet states with ft = A ± 1/2. Such dou- 
blets play a crucial role in explaining features of de- 
formed nuclei, including superdeformation and identical 
bands [|,[ll|. 

The symmetry of class II, referred to as spin symme- 
try, occurs when the difference of the scalar and vector 
potentials is a constant, Vs(p,z) — Vv{p,z) = S . The 
symmetry group is again SU(2) and its generators 0] 



S t = 
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(0) 



commute with the Dirac Hamiltonian. The Dirac eigen- 
functions in the spin limit satisfy 



ft = ±1/2 



(7) 



and form degenerate SU(2) doublets. Their wave func- 
tions are of the form Q 



*o 1 =a + i/ 2 : {5,0,/,/-} , 



(8b) 



where A = ft — p > is the eigenvalue of J z — S z . 
The upper components of the two states in Eq. ((SJ) form 
the usual non-relativistic spin doublet with a common 
radial wave function, g, an orbital angular momentum 
projection, A, and two spin orientations ft = A±l/2. The 
relativistic spin symmetry has been shown to be relevant 
to the structure of heavy-light quark mesons [|| . 

The Dirac Hamiltonian has additional symmetries 
when the scalar and vector potentials depend on differ- 
ent variables. The symmetry of class III occurs when the 
potentials are of the form Vs — Vs(z) and Vy = Vv(p)- 
In this case, the Dirac Hamiltonian commutes with the 
following Hermitian operator 

R z = [M + V s (z)]P± 3 +~, 5 p z , (9) 
where Ej = ( °£ ° \ The Dirac eigcnfunctions satisfy 



(10) 



A separation of variables is possible by choosing the Dirac 
wave function in the form 



where u, = Ui{p), h± = h±(z) and, for simplicity, we 
have omitted the label e from these wave functions. The 
Dirac equation then reduces to a set of two coupled first- 
order ordinary differential equations in the variable p, 

[d/dp-n/p} Ul (p)-[E-Vv(p)+e\u 2 {p) = (12a) 
[d/dp + n/p]u 2 {p) + [E-Vv{p)-e]u 1 (p) = (12b) 

and a separate set in the variable z 

[M + Vs(z) + d/dz] h 2 (z) = eh x {z) (13a) 
[M + V s (z)-d/dz]h 1 (z) = eh 2 (z) (13b) 

where h±(z) = h 2 (z)±hi(z). The separation constant, e, 
plays the role of a mass for the transverse motion and is 
determined from imposed boundary conditions. A special 
case within the symmetry class III, with Vs(z) — and 
e = ±y M 2 + pi, was considered for electron channeling 
in crystals For Vs(z) — 0, R z of Eq. ([9|), reduces to 
the transverse polarization operator relevant to studies 
of synchrotron radiation in storage rings and QED pro- 
cesses in magnetic flux tubes {e.g., e + e~ production and 
Bremsstrahlung) [l3j |. 

The symmetry of class IV occurs when the potentials 
are of the form Vs = Vs(p) and Vy = Vy(z). In this 
case, the following Hermitian operator 

R p = [M + Vs(j>)]ta-i0'ft(±Xp)3 (14) 

commutes with the Dirac Hamiltonian and the Dirac 
eigenf unctions satisfy 



(15) 



Again, a separation of variables is possible with the choice 
of wave function, 



a : itiw-, -bw + }/y/p , (16) 



where ^ = ^i(p) and w± = w±(z). The Dirac equation 
then reduces to a set of ordinary differential equations in 
the variable p, 

[d/dp - ft/p] ^(p) - [e + M + Vs(p) } UP) = (17a) 
[d/dp + n/p]Up) + ^-M-V s (p)]Up) = (17b) 

and a separate set in the variable z 

[E — Vv(z) — id/dz] w 2 (z) = ewi(z) (18a) 
[E -V v (z) + id/dz]w 1 (z) = ew 2 (z) (18b) 

where w±(z) — w 2 (z) ±wi(z). The quantum number, e, 
plays the role of an energy for the transverse motion. 
A particular selection of potentials within the symmetry 
class IV was encountered in the study of the Schwinger 
mechanism for particle-production in a strong confined 
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Vtfj : {u\h + , u 2 h-, Uih-, —u 2 h + } j yfp , (11) field (Vv(z) = ayz) [3, 15|. qq pair-creation in a flux 
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TABLE I: Conserved, anticommuting operators for Dirac Hamiltonians (H) exhibiting a super symmetric structure. 
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R 
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= f(H) 


V s (p, z) + V v (p, z) = A 


k © 


2(M + A - H)L 


{M + 


A - H) 2 


Vs(p,z) - V v (p,z) = H 


s z © 


2( M + H + H) 4 


(M + 




V s = Vs(z), V v = s f 


R z © 


(3± 3 {iJ zl5 [H - t 3 R\] - sx-(± ■ p) R z } 


J 2 (H 2 ~ 


Ri)+a 2 v Rl 


Vs = Sf,Vv = V v (z) 


r p ci 


£ 3 {i.L l5 [M - ± 3 R P ] -2*(±> p)0R P } 


Jz(Rp - 


- M 2 ) + a 2 s R 2 p 



tube (V s (p) = 0, e = ±^/M T +J?) [l6|, and the canon- 
ical quantization in cylindrical geometry of a free Dirac 
field (V s (p) = Vv(z) = Q,E 2 = M 2 + k 2 +p 2 z ) 0. 

Dirac Hamiltonians with selected external fields are 
known to be supersymmetric [!, H, H, EH . It is, therefore, 
natural to inquire whether a supersymmetric structure 
can develop within each of the above symmetry classes. 
The essential ingredients of supersymmetric quantum 
mechanics (lij are the supersymmetric Hamiltonian, Ti, 
and charges Q+, Q_ = Q+, which generate the super- 
symmetry (SUSY) algebra [H,Q±] = {Q±,Q±} = 0, 
{Q-iQ+} = Accompanying this set is an Hermitian 
Z 2 -grading operator satisfying [H,V] — {Q±,V} = 
and V 2 = 11. The +1 and —1 eigenspaces of V define the 
"positive-parity", H + , and "negative-parity", if_, sec- 
tors of the spectrum, with eigenvectors vI/W and 
respectively. The SUSY algebra imply that if <I>( + ) is an 
eigenstate of Tt, then also <J/(-) = is an eigen- 

state of Tt with the same energy, unless Q^^^ vanishes 
or produces an unphysical state, (e.g., non-normalizable). 
The resulting spectrum consists of pairwise degenerate 
levels with a non-degenerate single state (the ground 
state) in one sector when the supersymmetry is exact. 
If all states are pairwise degenerate, the supersymmetry 
is said to be broken. Typical spectra for good and bro- 
ken SUSY are shown in Fig. 1. Degenerate doublets, 
signaling a supersymmetric structure, can emerge in a 
quantum system with a Hamiltonian H, from the exis- 
tence of two Hermitian, conserved and anticommuting 
operators, R and B 



(19) 



[H, R] = [H, B] = {R, B} = 



Q. 

► 

H + H_ H + H_ 

good SUSY broken SUSY 

FIG. 1: Typical spectra of good and broken SUSY. The op- 
erators Q- and Q+ connect degenerate states in the H+ and 
H- sectors. 



The operator R has non-zero eigenvalues, r, which come 
in pairs of opposite signs. B 2 = B^B = f(H), is 
a function of the Hamiltonian. A ^-grading opera- 
tor, V r = R/\r\, and Hermitian supercharges Q\ = B, 
Q2 = iQ{P r can now be constructed. The triad of oper- 
ators Q± = (Qi ± iQ 2 )/2 and U = Q\ = f(H) form the 
standard SUSY algebra. In the present analysis, f(H) is 
a quadratic function of the Dirac Hamiltonian, H , and 
the relevant R and B operators are listed in Table I. 

In the pseudospin symmetry limit, the relevant opera- 
tor B, connects the doublet states of Eq. ([5]). The spec- 
trum, for each A ^ 0, consists of twin towers of pairwise 
degenerate pseudospin doublet states, with Q,\ = A — 1/2 
and f^2 = A + 1/2, and an additional non-degenerate 
nodeless state at the bottom of the Sli = A — 1/2 tower. 
Such nodeless states correspond in the non-relativistic 
nuclear deformed shell- model to the "intruder" states, 
[N, n 3 , A = N — n 3 ]ri = A + 1/2, which, empirically, 
are found not to be part of a doublet ll|. The latter 
property follows from the fact that a nodeless bound 
Dirac state satisfies the criteria of Eq. ([2]), hence has 
a wave function as in Eq. (|5a|) with g + , g, f ^ and 
f/g + > 0. Its pseudospin partner state has a wave 
function as in Eq. (|5bj) . The radial components sat- 
isfy Bg~ — [B — 2(A / p)f I g + ]g + , where B is defined in 
Eq. (|lb[) . This relation is satisfied, to a good approxima- 
tion, for mean-field potentials relevant to nuclei, and the 
r.h.s. is non-zero and, consequently, g~ 7^ 0. If so, then 
the partner state f|5b|) is also nodeless, but it cannot be a 
bound eigenstate since its radial components do not ful- 
fill the condition of Eq. @. Altogether, the ensemble of 
Dirac states with O2 — ^1 = 1 exhibits a supersymmetric 
pattern of good SUSY, as illustrated in Fig. (2a). 

In the spin symmetry limit, the relevant operator B 
connects the doublet states of Eq. |8]| . The spectrum, for 
each A ^ 0, consists of twin towers of pairwise degenerate 
spin-doublet states with f2i = A— 1/2 and Q2 = A + 1/2. 
None of these towers have a single non-degenerate state. 
This follows from the fact that, in view of Eq. @, a 
nodeless bound state has a wave function as in Eq. (|8a[) 
with g, f,f~^0 and g/f~ > 0. Its spin partner has 
a wave function as in Eq. (I8bj) . The radial components 
satisfy Af+ = [A - 2 (A/ p)g/ f~]f~ , where A is defined 
in Eq. (fTa|). For relevant potentials the r.h.s. of this 
relation can vanish, hence / + has a node. Therefore, the 
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[411] 1/2 [400] 1/2 [402] 3/2 (a) 

[431] 1/2 [420] 1/2 [422]3/2 [411]3/2 [413]5/2 [402]5/2 [404]7/2 

[440] 1/2 [431]3/2 [422]5/2 [413]7/2 [404]9/2 

A = A=l A = 2 A = 3 A = 4 A = 5 

[400] 1 /2 (b) 

[420] 1/2 [411]3/2 [411]l/2 [402]5/2 [402]3/2 

[440] 1/2 [431J3/2 [431J1/2 [422]5/2 [422]3/2 [413]7/2 [413]5/2 [404]9/2 [404]7/2 
A=0 A=l A=2 A=3 A=4 

FIG. 2: Grouping of deformed shell-model states [N = 4, n3,A]fi, exhibiting a pattern of (a) good SUSY, relevant to the 
pseudospin symmetry limit, and (b) broken SUSY, relevant to the spin symmetry limit. TV and 713 are harmonic oscillator 
quantum numbers. A (A) is the orbital (pseudo-orbital) angular momentum projection along the symmetry z-axis. 



spin-partner of a nodeless state is not nodeless and can 
be a bound eigenstate, since the restrictions of Eq. (jSJ) 
do not apply. Altogether, the ensemble of Dirac states 
with O2 — S^i = — 1 exhibits a supersymmetric pattern of 
broken SUSY, as illustrated in Fig. (2b). 

Within the symmetry class III, a supersymmetry is 
obtained for Vv(p) = o-v I P an d Vs{ z ) arbitrary. The 
energy eigenvalues are E^ n = \e\/y/l + tty/K + l) 2 
(rip = 0, 1, 2, . . .), with 7 = y/tl 2 - 0$. From Eqs. (JT3J) 
we see that if [hi(z), /i2(z)] are solutions with e > 0, then 
[hi(z), — Ji2(z)] are solutions with — e < 0. Accordingly, 
the doublet wave functions are as in Eq. (fTTj) . with the 
replacements, u, 1— ► uf\p) for VE^ 6 ' n , and Ui 1— > u\ e \p), 

h± h-> -h T (z) for ^ { ~%. For n p > 1, the states ^% 
are degenerate. For n p = only one state is an accept- 
able solution, which has e > (assuming ay < 0) and 
is annihilated by the relevant operator B. For each 
and e the spectrum resembles a supersymmetric pattern 
of good SUSY, with the towers H + {H- ) of Fig. 1 corre- 
sponding to states with e > (e < 0). 

Within the symmetry class IV, a supersymmetry is ob- 
tained for Vg(p) = ots/p (ens < 0) and Vy(z) arbitrary. 
The allowed values are e = ±Myl — a s /(n p + -y) 2 
(ri p = 0, 1, 2, . . .), where 7 = y/Q, 2 + a|. From Eqs. ^ 
we see that if [u>i(z), W2{z)} are solutions with e > 0, then 
[wi(z), — W2{z)\ are solutions with — e < and the same 
energy, E. Accordingly, the doublet wave-functions are 
as in Eq. (fTB|) . with the replacements, ft 1— > ft (fi) f° r 
< } ,a> and ft » ^(p), w ± » -w T (z) for &~%. For 

?i p > 1 the states are degenerate. For n p — only 

one state, with e > 0, is an acceptable solution, which is 
annihilated by the relevant operator B. Again, for each Q, 
and e the resulting spectrum resembles a supersymmetric 
pattern of good SUSY. 

In summary, we have considered classes of symmetries 
and related supersymmetries of Dirac Hamiltonians with 
cylindrically-deformed scalar and vector potentials. The 



symmetries arise when the potentials obey a constraint 
on their sum or difference, or when they depend on dif- 
ferent variables. The known pseudospin and spin sym- 
metry limits are by themselves supersymmetric. Addi- 
tional supersymmetries arise when one of the potentials 
has a 1/p dependence and the second potential depends 
on z. It is gratifying to note that some of the indicated 
(super)symmetries are manifested empirically, to a good 
approximation, in physical dynamical systems. 
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